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Abstract 

vQ ' An appropriate gauge-invariant framework for the treatment of magnetized curvature and en- 

^D ' tropy modes is developed. It is shown that large-scale magnetic fields, present after neutrino 



decoupling, affect curvature and entropy perturbations. The evolution of different magnetized 
modes is then studied across the matter-radiation transition both analytically and numerically. 

O I From the observation that, after equality (but before decoupling) the (scalar) Sachs- Wolfe con- 

tribution must be (predominantly) adiabatic, constraints on the magnetic power spectra are 

C^ ' deduced. The present results motivate the experimental analysis of more general initial con- 

ditions of CMB anisotropics (i.e. mixtures of magnetized adiabatic and isocurvature modes 
k> , during the pre-decoupling phase). The role of the possible correlations between the different 

Vh ' components of the fluctuations is partially discussed. 
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1 Introduction 

The observational evidence inferred from the spatial variations of the Cosmic Microwave Back- 
ground (CMB) temperature suggests that the curvature perturbations that are present prior to 
decoupling (but after equality) must be predominantly adiabatic and Gaussian with an (almost) 
scale-invariant spectrum [1, 2, 3]. A mode is defined to be (predominantly) adiabatic if the fluctua- 
tions of the entropy density are subleading after matter-radiation equality (but prior to decoupling) . 
An exactly adiabatic mode of curvature fluctuations may be produced, for instance, by single-field 
inflationary models. If this is the case, the scalar fluctuations of the inflaton field (with wavelength 
larger than the Hubble radius) are transferred to curvature fluctuations whose value stays (almost) 
frozen from the early stages of the evolution of the Universe until the pre-decoupling epoch when the 
minute spatial variations of the geometry are imprinted, via the Sachs- Wolfe effect, on large-scale 
temperature inhomogeneities. 

In addition to adiabatic fluctuations there may be one or more non-adiabatic modes sometimes 
named (with a somehow confusing terminology) isocurvature modes. The specific number of (reg- 
ular) isocurvature modes (up to four) is determined by the physical properties of the plasma where 
curvature fluctuations evolve. A mode, in the present terminology, simply means a consistent so- 
lution of the governing equations of the metric and plasma fluctuations, i.e. a consistent solution 
of the perturbed Einstein equations and of the lower multipoles of the Boltzmann hierarchy. 

Consider as an example the physical system formed by two fluids, i.e. radiation and pressureless 
matter. This is a pretty good approximation in the description of pre-decoupling physics if the 
radiation is identified with the sum of photons and neutrinos and the pressureless fluid accounts 
for the cold dark matter (CDM) component. If this is the case, the plasma suports not only an 
adiabatic mode, but also a non-adiabatic mode for which the entropy fluctuations of the CDM- 
radiation system do not vanish [4, 5]. It is therefore possible that the pre-decoupling curvature 
fluctuations have a dominant (adiabatic) contribution and a subdominant (entropic) correction. 
Sticking to the simplified situation described above, the large-scale contribution to the angular 
power spectrum would be determined by three pre-decoupling power spectra, i.e. the spectrum of 
the adiabatic mode, the spectrum of the entropic mode and the spectrum of their cross-correlation. 

This possibility leads to a rich parameter space and to an interesting phenomenological frame- 
work (see, for instance, [6, 7, 8]). As an example in Ref. [8] a recent likelihood analysis determined 
the isocurvature fraction allowed by WMAP data (as well as by large scale structure data) by 
taking into account both the adiabatic and the entropic contribution of the CDM-radiation system. 

The study of isocurvature and entropy modes is often conducted within a two-fluid system 
because of the inherent mathematical complications related with the possible proliferation of non- 
adiabatic modes. It will be demonstrated here that the presence of a large-scale magnetic fields 



during the pre-decoupling epoch affects the whole system of curvature and entropy perturbations 
both quahtatively and quantitatively. 

If the origin of large-scale magnetic fields is primordial (as opposed to astrophysical) it is 
plausible to expect that the pre-decoupling plasma is effectively magnetized. Different studies have 
been devoted, in recent years, to constrain a hypothetical primeval magnetic field from CMB studies 
(see [9] for a short review and [10] for a more extended discussion). Indeed, large-scale magnetic 
fields that are present in the pre-decoupling epoch may affect the polarization power spectra (both 
the EE, BB and EB angular power spectra) either directly [11, 12, 13] or through the Faraday 
effect [14, 15, 16, 17, 18, 19] that may rotate an initial E mode (produced thanks to the polarized 
nature of Thompson scattering) into a B mode. Barrow and Subramanian, in a series of papers 
[20, 21, 22, 23], analyzed the vector modes (and the polarization power spectra) induced by large- 
scale magnetic fields. Recently [12] a numerical approach has been developed to analyze vector and 
tensor modes. Vector, as well as tensor modes have been discussed also in [11]. While the effect of 
vector modes will be crucial at small scales (i.e. i > 1000) the contribution of tensor modes may 
be relevant also for larger scales (i.e. £ < 100). 

Large-scale magnetic fields induce vector and tensor modes but also various scalar modes whose 
classification and analysis will be one of the themes of the present investigation. More specifically, 
large-scale magnetic fields affect the usual adiabatic and entropic modes [24, 25] that may be present 
in the pre-decoupling epoch. In fact, the divergence of the magnetohydrodynamical Lorentz force 
influences the tight-coupling evolution of photons and baryons. Moreover, the magnetic fields 
gravitate. Therefore, they contribute to the Hamiltonian constraint and (in a weaker fashion) to 
the momentum constraint. In [26], an interesting covariant approach for the analysis of cosmological 
perturbations induced by large-scale magnetic fields has been proposed and further developed in 
[27]. The approach used in [27] was to assume a background magnetic field which uniform but, 
at the same time, sufficiently weak not to break the isotropy of the background geometry. The 
present approach is more similar to the one of Ref. [28] where the magnetic field has been taken 
fully inhomogeneous. Moreover, the problems related to the present analysis concern the specific 
modifications induced by large-scale magnetic fields on the adiabatic and non-adiabatic modes that 
are normally used to set initial conditions for CMB anisotropics. This goal requires a full description 
of the magnetized photon-baryon system in the presence of neutrinos and CDM particles both before 
and after decoupling. The strategy would then be to study how large-scale magnetic fields affect 
the system of equations describing adiabatic and non-adiabatic modes before the time of radiation- 
matter equality, i.e. deep in the radiation-dominated epoch. The magnetized solutions will be then 
evolved both analytically and numerically across the radiation-matter transition and up to the time 
of photon decoupling. 

As mentioned, large-scale magnetic fields are assumed here to be fully inhomogeneous and char- 



acterized by an appropriate two-point function. In this way the spatial isotropy of the background 
geometryis not broken and the Hne element can then be written, in the spatially flat case, as 

ds'^ = gf.udx'^dx'' = a^{T)[dT^ - dx'^], (1.1) 

where r is the confornial time coordinate and a{T) is the scale factor. Defining the appropriately 
rescaled magnetic fields as B = a^B (see section2), the relevant two-point function can be written 
(see the appendix for a summary on the conventions) as 

{B,{k,T)B\p,T)) = ^P^(k)PBik,T)6^^Hk+p), (1.2) 

where PB^k) is the magnetic field power spectrum and P- (k) = {6j —kik^/k"^). Large-scale magnetic 
fields will be treated here in a magnetohydro dynamical (MHD) approximation [29] in curved space- 
time. This description is justified by different considerations. In fact, before decoupling, not only 
the plasma is globally neutral but the value of the Debye scale much smaller than the Hubble 
radius and it is also much smaller than the electron and photon mean free path [10]. As far as 
the magnetic properties of the system are concerned, the plasma behaves, in this situation, as a 
single (intrinsically charged but globally neutral) fiuid where the total Ohmic current is solenoidal 
and the Ohmic electric fields are suppressed by the finite value of the conductivity. The same 
kind of MHD description has been used also in the analysis of the vector and tensor modes of the 
geometry (see, for instance, [11, 12] and [20, 21]). It should be borne in mind, however, that it 
would be incorrect to use the same kind of MHD approximation in the analysis of Faraday effect 
(as sometimes assumed). Faraday effect concerns the propagation of electromagnetic waves in a 
cold (or warm) plasma and cannot be treated within MHD which is valid for typical frequencies 
much smaller than the plasma frequency and typical scales much larger than the Debye scale. So 
in full analogy with the discussion of Faraday effect in the laboratory (see for instance [30, 31]) a 
two-fluid description is mandatory [17]. 

A direct consequence of the MHD approach is that magnetic flux and magnetic helicity will 
be, to a good approximation, conserved in the limit of large conductivity [29]. In the Friedmann- 
Robertson- Walker background (1.1), flux conservation implies that the magnetic field strength 
increases as a~^ as we go back in time. At the present time magnetic fields are observed, at the 
ji G level, in galaxies [32], clusters [33, 34] and in some superclusters [35]. Even if the fields have 
different properties in diverse classes of structures and even if the determination of magnetic fields 
in superclusters is still in progress, it is rather plausible to admit that we live in a sort of magnetized 
Universe [36, 37]. In the case of galaxies, if the magnetic field present prior to gravitational collapse 
was of the order of 0.1 nG, compressional amplification alone could have produced the observed 
magnetic field even without the need of a dynamo action. In recent years interesting progresses 
have been made in the analysis of the dynamo problem especially in connection with the role of 
global conservation laws of magnetized plasmas [38, 39]. Therefore, the protogalactic field could 



be much smaller than 0.1 nG and still be explain the observed magnetic fields in galaxies, while in 
this framework cluster magnetic fields are related to the way dynamo saturates. Aside from this 
consideration, if we blue-shift a protogalactic magnetic field (coherent over a typical Mpc scale) of 
0.1 n G at the epoch of photon decoupling we will end up with a magnetic field of mG strength 
that can potentially affect CMB anisotropics. The mechanisms that may produce magnetic fields 
will not be specifically examined and we refer the reader to a recent revies [37]. However, most 
of the mechanisms predict a stochastically distributed magnetic field that does not break spatial 
isotropy. It will then be interesting to scrutinize the implications of a pre-decoupling field on the 
scalar CMB anisotropies. In this sense the ambition of the present approach is a model-indepedent 
approach that is similar to recent attempts of constraining mixtures of adiabatic and non-adiabatic 
modes in spite of their specific origin (see [6, 7, 8] and references therein). 

To be more specific, the strategy will be to treat the magnetic fields as a further element of 
the primordial plasma that may contain both adiabatic (or entropic) and non-adiabatic modes. 
In doing a gauge invariant description will be adopted since, as it will be illustrated, a simpler 
description of the matter-radiation transition can be achieved. To say anything concerning CMB 
anisotropies we have to estimate what is the value of curvature perturbations after equality (but 
before decoupling) and in the long wavelength limit ^. Such an estimate, is, unfortunately, not 
available in the present literature. It is then appropriate to introduce the description developed 
in this paper by defining the density contrast on uniform curvature hypersurfaces, conventionally 
denoted by C [5, 40, 41] 

C = -H^, Pt = EPi' p[ = -m{pt+pt), (1.3) 

Pt i 

where the prime denotes a derivation with respect to r and 7i = (In a)'; pt is the (total) background 
fiuid density while pt is the pressure. Similarly, in Eq. (1.3), 6pt = 5p + 5pB, Sp is the fiuid density 
fiuctuation in the uniform curvature gauge [42, 43] and 5pB is the fully inhomogeneous magnetic 
energy density which is gauge-invariant and quadratic in the magnetic field intensities. If the fiuid 
has different consituents (for instance, in one of the simplest cases, dusty matter and radiation), 
then the total density contrast on uniform curvature hypersurfaces, i.e. 6pt can be written as the 
weighted sum of the individual fiuid perturbations and of the magnetic field density contrast, i.e. 

C = e4o + Cb, = -H^, Cb = -H^. (1.4) 

Y Pt Pi Pt 



•^In the present paper we will often refer to the "long wavelength limit" of some expressions or equation. By this 
we simply mean that the relevant wavelengths are the ones whose typical size is larger than the Hubble radius at the 
corresponding epoch. If we define the comoving wave-number k, the long wavelengths, in the present terminology, 
are the ones satisfying fcr < 1. 



Given the form of Eqs. (1.4) the relative fluctuations in the entropy density ? can be defined as 

5ij = ^ = -3(Ci - Cj), cSiB = ^ = -3(Ci - Cb). (1.5) 

In this language, the adiabatic mode corresponds to the situation where Cij = and CiS = for each 
fluid of the mixture. The other important variable treated here will be gauge-invariant curvature 
fluctuation ^ usually denoted by TZ, and related to ( as 

12 7rGa^(pt + pt) 

where ^ is the (gauge-invariant) Bardeen potential. In the long wavelength limit, the second term 
at right hand side of Eq. (1.6) can be neglected and C is approximately equal to TZ. When the 
plasma is constituted by a mixture of relativistic fluids, the evolution of C can be written, in the 
long wavelength limit , as 

C' = '^Spn^d + —^(c^-\)5pB, (1.7) 

Pt + Pt Pt + Pt\ 3/ 

where (5pnad is the non-adiabatic pressure density variation accounting for the contribution of en- 
tropy fluctuations (see Eq. 2.31). Equation (1.7) contains, at the right hand side, not only the 
contribution of 5pnad but also the effect of the inhomogeneous magnetic field. 

The logic followed in the present investigation will then be to provide a fully gauge-invariant 
system of equations where the evolution of magnetized curvature and entropy modes can be char- 
acterized as regular solutions deep in the radiation-dominated regime. Then the various modes will 
be followed across the radiation-transition both analytically and numerically. 

In Section 2 a gauge-invariant system of governing equations will be proposed for the simulta- 
neous description of magnetized adiabatic and entropy modes. In Section 3 the analysis of various 
magnetized modes will be presented in the realistic situation where the plasma is constituted by a 
mixture of baryons, neutrinos, CDM and photons. The fate of the regular magnetized modes across 
the radiation-matter equality will be discussed in Section 4. The precise amplitude of the curvature 
perturbations before decoupling will then be computed for the wavelengths that are larger than 
the Hubble radius after equality. In Section 5 some more phenomenological considerations on the 
properties of the Sachs- Wolfe plateau will be discussed with the aim of providing a possible starting 
point for more specific numerical discussions. Finally, Section 6 contains the concluding discussion 
and a summary of the possible developments of the present investigation. A series of technical 
results have been collected in the Appendix in connection with the proper parametrization of the 
various power spectra arising in the problem. 

^In the present paper curvature fluctuations are identifled with the fluctuations of the spatial curvature in the 
comoving orthogonal gauge [44, 45]. The gauge-invariant generalization of this quantity is usually denoted by 7?.. 
Clearly, in different gauges, the physical meaning of 72. may be different. 



2 Magnetized curvature fluctuations 

2.1 General considerations 

The system of magnetized curvature perturbations can be obtained by considering the case of 
a gravitating charged fluid which is, however, globally neutral for scales larger than the Debye 
length. Around the time of decoupling, for a typical temperature of 0.3 eV and for an electron 
density ne = 10^ cm~^ the Debye length is of the order of 10 cm, i.e. much smaller than the Hubble 
radius at the corresponding epoch. The background equations (and the background metric (1.1)) 
are not affected by the presence of magnetic fields and have the standard form in terms of the 
(total) energy and pressure densities, i.e. pt and pt- 

ji^ = ^a^p^^ n''-n' = 4TTGa\pt+pt), p[ + 3nipt+Pt) = 0. (2.1) 

Denoting with 6s the first-order (scalar) fluctuation of the corresponding quantity, the perturbed 
Einstein equations can be written as 

4a;: = 8^G(4r; + r;), (2.2) 

where Q^ is the Einstein tensor while (JgT^ and 7^ are, respectively, the energy-momentum tensors 
of the fluid perturbations and of the electromagnetic sources. The perturbed covariant conservation 
equations and by the equations for the Maxwell fleld strength can instead be written as 

<5s(V^T,^)+V^T,^ = (2.3) 

V^F^^ = inf. (2.4) 

In Eqs. (2.2)-(2.4), V^ denotes the covariant derivative defined from the metric g^u defined in Eq. 
(1.1). In Eq. (2.4) f = (TcF^^'^u^ is the Ohmic current. The (total) velocity field, w^, obeys, as 
usual g^^'^u^Uy = 1. Recalling that Fqi = aP'Si and Fij = —a^eij^B^ the evolution equations for the 
rescaled electric and magnetic fields become: 

V-^ = 0, V-i? = 0, (2.5) 

dB - - 

'd- = -VxE, (2.6) 

- - ^ dE 

VxS = 47rJ + ^— , (2.7) 

OT 

J = a{E + vxB), (2.8) 

where 

E = c?e, B = c?B, J = a^j. (2.9) 



The one-fluid MHD description adopted here [29] amounts then to the requirement that both the 
Maxwell fields and the induced Ohmic current are solenoidal 

V.^ = 0, V-i? = 0, V-J = 0, (2.10) 

which implies that 

J=^ , E=— vxB. (2.11) 

Since we will be interested in effects occurring over length-scales much larger than the Debye scale, 
the displacement current can be safely neglected in Eq. (2.8). Note that v appearing in Eq. (2.11) 
is the bulk velocity of the plasma. In MHD, the bulk velocity of the plasma is essentially the center 
of mass velocity of the baryon-electron system. Therefore, v coincides, in practice with the baryon 
velocity field. In Eq. (2.11), a is the conductivity. In a relativistic plasma a ~ T/acm while in 
a non-relativistic plasma a ~ (T/acm){T/me)^''^ where rUe is the electron mass and T is common 
temperature of electrons and baryons. 

The evolution equations for the coupled system defined by Eqs. (2.2), (2.3) and (2.4) can then 
be derived. For the present purposes it is useful to employ the uniform curvature gauge [42, 43] 
where the perturbed elements of the metric tensor are 

Ssgoo = 2a^(j), dsQio = -a^diA. (2.12) 

The covariant conservation of the total energy-momentum tensor, i.e. Eq. (2.4) implies then two 
separate equations one for the total velocity field and the other for the density fiuctuation. These 
two equations can be written as 

5sp' + (Pt + Pt)e + m{5,p + M = ^, (2.13) 

{6 + V^A)' + n{l - 3c2)(0 + V^^) + ^^ - ^^^ + VV = ^:y''-^l (2.14) 

Pt + Pt Pt + Pt a^{pt + Pt) 

where 6 = div'^ is the divergence of the (total) velocity field and where Cg = p[/ p[ is the (total) 
sound speed. To derive Eqs. (2.13) and (2.14) from Eq. (2.3), the perturbed components of the 
energy-momentum tensor of the fluid have been parametrized as 

5,T^ = 6sp, 5Ti = -SsP + n|, (2.15) 

6,T^ = (pt + Pt)v\ 6sTJ' = -{pt+ pt){v^ + ^^A). (2.16) 

In Eq. (2.16) the term 11^ represents the contribution of the anisotropic stress. The energy- 
momentum tensor of the electromagnetic sources is, instead, 

V = ^ {-F,c.F^^ + -/,F^pF^P\ . (2.17) 



Equation (2.2) implies two constraints (stemming from the (00) and (Oi) components). The Hamil- 
tonian and momentum constraints are then, in the ofF-diagonal gauge, 



v'^[n<t) + {n^ - n')A] = -^T:Ga\pt + pt)e + 5VbI 



where 



5PB 



B' 



5Vb 



V -{E X B) 



(2.18) 
(2.19) 

(2.20) 



are, respectively, the magnetic energy density and the Poynting vector. Notice that 6Vb, in the 
one-fluid MHD approximation, contains two terms: the first is suppressed by the conductivity, the 
second is of higher-order. In fact, using the second relation of Eq. (2.11), 



SVt> 



V-{J X B) V -[{vx B)x B] 



The (ij) component of Eq. (2.2) can be written as 

{n^ + 2n')(i) + He/)' + ^v^{A' + 2nA + 



(2.21) 



5' 



-did^[A' + 2nA + . 



where 



A7rGa\6sP + 5pB)Si-U{-Ul], 



n^ 



SPB = — , 



47ra4 



B,B^ - ^B^Sj 



(2.22) 



(2.23) 



The electromagnetic contribution to the Hamiltonian constraint (2.18) contains also a purely 
electric term. However, by virtue of Eq. (2.11) the electric contribution is suppressed with the 
respect to the magnetic one. Furthermore, terms containing the product of a quadratic (magnetic) 
contribution and of a metric (or velocity) fluctuation are of higher order and will then be neglected. 



2.2 Gauge-invariant governing equations 

The system formed by Eqs. (2.13), (2.14) and (2.22) (supplemented by the constraints (2.18) and 
(2.19)) mixes the four variables 0, A, 5sP and 6. The gauge-invariant generalization of the total 
density contrast on uniform curvature hypersurfaces, i.e. 6sP is given by C, i.e. 



c 



^sPt^ 

— r^' 
Pt 



(2.24) 



where 5sPt = SsP + SpB ■ On the other hand, the gauge-invriant generlization of (f) is related with 
the curvature fluctuations on comoving orthogonal hypersurfaces, denoted by TZ: 



n 



n^-n 



i<t^- 



(2.25) 



Finally, the variable A is related with one of the two Bardeen potentials, usually denoted by ^: 

qt = -nA, (2.26) 

The variables (", ^ and ^ are all gauge-invariant and while they have been introduced, for simplicity, 
in the uniform curvature gauge, they have interesting physical interpretations also in different 
gauges. For instance, ( can be interpreted as the curvature perturbation on constant density 
hypersurfaces. This interpretation is clear by looking at the gauge-invariant form of the Hamiltonian 
constraint. Inserting Eqs. (2.24), (2.25) and (2.26) into Eqs. (2.18) and (2.19) the gauge-invariant 
forms of the Hamiltonian and momentum constraints can be obtained: 

UirGa-'lpt +pt) 

\/^(n + ^)=ne + <5sVb. (2.28) 

pt + pt 

In the gauge-invariant approach presented here, the Hamiltonian constraint (2.27) fixes the relation 
between ^ and TZ. If we are interested in the evolution of the system for length-scales much 
larger than the Hubble radius we will have that the term containing V^^ is sub-dominant and 
the constraint is satisfied if C — ^- The momentum constraint (2.28) can be used to compute the 
divergence of the total velocity field 9 once TZ and ^ are known. Using Eqs. (2.24), (2.25) and 
(2.26) into Eq. (2.13), and recalling Eqs. (2.1), the following (gauge-invariant) evolution equation 
for (" can then be obtained 



c " 



SPnad - ( Cs - - ) (JPE 



I + Ell (2 29) 

Pt + Ptl^ '"-^ V ^ 3/ ■ "J 3 3a4(pt + pt) ' ^ ■ ^ 

In Eq. (2.29), the pressure fluctuation of the fluid component has been split as 6sP = c^<JsP + 5pnad; 
where (5pnad is the term proportional to the intrinsic entropy fluctuations of the fluid and Cg = p[/p[ 
is the total sound speed. This decomposition is usual when dealing with mixtures of relativistic 
fluids even in the absence of magnetic fields. The explicit form of <5pnad can be evaluated once the 
chemical composition of the fluid is specified. By introducing the individual sound speeds of the 
various fluids ^ and the relative fluctuations of the specific entropy in a generic pair of fluids, i.e. 

Csi = 4 5ij = -3(Ci-Cj), (2.30) 

Pi 

the general form of 5pnad can be written as 

^Pnad = -^^ Y^ PiPj(Csi - Csj)'5ij. (2.31) 

The indices appearing in Eq. (2.31) (not to be confused with covariant spatial indices) label the 
different components of the fluid and, by definition, Sij = — 5j i (see Eq. (2.30). The gauge-invariant 
^Clearly, by definition of total sound speed c^ = X].(Pi/pt)Csi- 



density contrasts of each fluid, i.e. Ci) can be evaluated, for instance, in the uniform curvature gauge 
where they correspond to 

Ci = -W^, Pt = Y.P^ (2-32) 

Pi I 

as already introduced in Eq. (1.4). The variable C, that appears in the Hamiltonian constraint 

(2.27) can be also written in terms of the individual C,\ supplemented by tha magnetic contribution 

Cb: 

C = E4Ci + CB, Cb = -H^. (2.33) 

Y^t Pt 

In connection with Eq. (2.29) the following remarks are in order: 

• even if the intrinsic (non-adiabatic) contribution vanishes, the evolution of C, may change over 
scales larger than the Hubble radius since its evolution is directly sourced by the magnetic 
energy density; 

• this new (entropic) contribution vanishes when the Universe is dominated by radiation but 
becomes relevant as soon as the speed of sound Cg 7^ 1/3 so that, for instance, it is relevant 
in the context of the radiation-matter transition; 

• by virtue of the momentum constraint (2.28) and of the Hamiltonian constraint (2.27), the 
contribution containing 9 at the left hand side of Eq. (2.29) is proportional, to leading order, 
to V^C &iid to V^^ (both contributions are then small on super-Hubble scales); 

• as previously remarked in deriving Eq. (2.28) the term E- J appearing in Eq. (2.29) is strictly 
zero in the superconducting limit and negligible for finite conduvctivity. 

The last statement can be easily verified since, in the MHD approximation, the relations of Eq. 
(2.11) imply the following chain of equalities 

E^./=-^''-^^^-(„-xg)^/='^'-^'''-'°^^''-'^^^''. (2.34) 

a Ana An 

Now, the first term in the second equality of Eq. (2.34) vanishes in the limit of large conductivities 
or, more precisely, in the limit a/H S> 1 where H = a7i is the Hubble parameter. The second 
term in the last equality of Eq. (2.34) is of higher order. To complete the discussion we have now 
to study the relations implied by Eq. (2.22) in terms of (", IZ and ^ . Using Eqs. (2.24), (2.25) and 
(2.26) into Eq. (2.22) we can extract the trace-free (gauge-invariant) contribution and apply the 
differential operator did^: 



V 



'^ + * + ^r^(^ + * 

3(/)t+Pt)vW 



-l^^[v^n + v^n]. (2.35) 

3(pt +Pt) 



10 



On the other hand, the trace-full contribution stemming from Eq. (2.22) simply leads to the 
gauge-invariant evolution of TZ, i.e. 



n' « 



Pt + Pt 



Spns.d - ( Cs - 2 ) (5pB 



''^^'^ +1^^{V'U + V'U). (2.36) 



47rGa2(pt + pt) 3pt + pt 



Equation (2.29) has the same physical content of Eq. (2.36). The reason is that for typical 
wavelengths larger than the Hubble radius, Eq. (2.27) implies that Q c^TZ. This statement can be 
verified by showing that, for instance, Eq. (2.36) reduces to Eq. (2.29) when the constraints (2.27) 
and (2.28) are taken into account. In fact, using Eq. (2.27) to express TZ in terms of C, and taking 
the first time derivative we can substitute TZ' in Eq. (2.36). This manipulation will give rise to a 
bunch of terms: the first one will be essentially C,' while the other contain the laplacian of ^ and 
its derivatives. Then we can use Eqs. (2.28) and (2.35) to get rid of ^'. The resulting equation 
will exactly be Eq. (2.29) which was obtained from the covariant conservation equation. 

Equations (2.29) and (2.36) suggest that the magnetic field acts, effectively, as a non-adiabatic 
variation of the pressure density. Therefore it is useful, sometimes, to think that 

<5Pnad = 5pnad " (Sc^ - 1)Cb, Cb = „/^_f , , (2.37) 

3(pt +Pt) 

where Cb now obeys 

C^+H(1-3c2)Cb = 0. (2.38) 



3 Magnetized curvature perturbations before equality 

3.1 Physical scales 

In the absence of magnetic fields, the commonly employed strategy is to use Eqs. (2.29) and (2.36) 
for the estimate of the values of curvature perturbations in the long wavelength limit after the 
equality (denoted, in what follows, by Tcq) but before decoupling (denoted by Tdec)- The magnetic 
contribution to Eqs. (2.29) and (2.36) change the evolution of C, and TZ, as discussed in the previous 
section. Moreover, the contribution to the anisotropic stress arises both from the collisionless 
neutrinos and from the magnetic fields. In the absence of collisionless neutrinos the magnetic fields 
can only induce, during radiation, a singular mode since the magnetic anisotropic stress cannot 
be counterbalanced by the anisotropic stress of the fluid. The solutions of the gauge-invariant 
system of section 2 will now be studied in the limit r <C Teq. Then the different magnetized modes 
(both adiabatic and non-adiabatic) will be followed through the matter radiation transition with 
the purpose of computing in detail the value of the curvature perturbations before decoupling and 
for wavelengths larger than the Hubble radius at the corresponding epoch. The radiation-matter 
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transition will be modeled by the interpolating scale factor 

(3.1) 



a(T 



-%2- 



which is a solution of Eqs. (2.1) when the plasma contains both pressureless matter and radiation. 
In Eq. (3.1), Oeq is the value of the scale factor at equality. The values of aeq and ri are fixed by 
the critical fractions of matter and radiation. In particular ^, 

Oen Tir 4.15x10"^ 2 fo^ f h'^Tl^,\-^ , , 

_zS. = _ ~ J-, = _^ ~ 287 48 — f3 2) 

ao ^n /i^S^m ' i/oV^m ' V 0.134/ ■ ^ ' ' 

and r^m and 0.^ are the present values of the critical fractions of matter and radiation, 

r^jn ^ r^c + r^b, rir = o.^ + o.^,. (3.3) 

The redshift of decoupling will be taken to be 1 + Zdec = oJec — HOO; implying from Eq. (3.2) that 
Tdec — 2.36 Teq (where Teq = (v2 " l)''"i)- The numerical values of the matter and radiation density 
fractions will be taken to be 

/i^Hin ~ 0.134, h^Tt, ~ 4.15 X 10"^ (3.4) 

with 

/i^Hc ~ 0.111, h^Vt^ ~ 0.023, h^Tl^ ~ 2.47 x 10"^ h^Vty ~ 1.68 x 10"^ (3.5) 

Concerning the set of parameters adopted in Eqs. (3.4) and (3.5) few comments are in order. The 
recent WMAP three year data [3], when combined with "gold" sample of SNIa [48] give values 
compatible with the ones adopted here. However, if WMAP data alone are used we would have, 
instead, a slightly smaller value for the critical fraction of matter, i.e. h?VL„i ~ 1.268. The values 
of the parameters chosen here are just illustrative in order to model, in a realistic fashion, the 
radiation-matter transition. Therefore, any difference in h?Q.^i will just affect the absolute value of 
Ti (and, consequently, of Teq) but will not change the features of the forthcoming numerical analysis 
whose purpose is just to corroborate the analytical results. A final comment on neutrino masses 
is appropriate. If neutrino masses are in the range of the meV (i.e. smaller than the temperature 
at ^dec)) they will be non-relativistic today. However, they will be relativistic around z^^c aiid will 
then contribute to the radiation fraction Or (and not to the matter fraction). In the following 
(when writing the evolution equations of the neutrino component) neutrino masses will then be 
neglected since our analysis applies for r < Tdcc- In setting the initial conditions (i.e. for r <C Teq) 
the neutrinos are effectively massless since the temperature of the radiation plasma is higher than 

-'dec 

^The scale factor will be normalized in such a way that ao — 1 and the factor h, i.e. the indetermination in the 
Hubble constant, will be taken 0.73 for numerical estimates. 
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3.2 Mixtures of four fluids 

For r <^ Tcq (but always after neutrino decoupling, i.e. T < 1 MeV) the radiation component 
dominates and, within the gauge-invariant approach presented in section 2, it is particularly sim- 
ple to set adiabatic or isocurvature initial conditions and to solve the equations deep inside the 
radiation-dominated phase. Few relevant cases will now be illustrated in view of the numerical 
analysis reported in section 3. When the plasma is formed by baryons, photons, CDM particles 
and massless neutrinos, the gauge-invariant density perturbations of each fluid obey 



3' '^ 3 



cL = -?, c; = --^, (3.6) 



C = -y, C = -f, (3.7) 

where the subscripts b, 7, c and ly refer, respectively, to the baryon, photon, CDM and neutrino 
components of the plasma. Since the divergences of the peculiar velocities carry a spatial gradient, 
they are small for modes that are larger than the Hubble radius prior to equality. Adiabatic initial 
conditions amount then to the choice that all the density contrasts are equal and constant, i.e. 

Cb = Cc = C7 = C;^ = constant. (3.8) 

Equation (3.8) implies that all the non adiabatic pressure density variation of the fluid vanishes 
exactly by virtue of Eq. (2.31). If the initial conditions for the gauge-invariant density contrasts are 
different from Eq. (3.8), then the corresponding solution is said to be non-adiabatic. In fact, any 
deviation from the initial conditions (3.8) necessarily leads to a non-vanishing entropy perturbation. 
For instance, if (^ 7^ C^, S^y ^ Q. This is the situation of the so-called neutrino isocurvature mode. 
Yet another example could be the one of the CDM-radiation mode where 5^c 7^ and Syc i^ ^■ 
The evolution equations of the velocity fields of the CDM and of the neutrinos are that appear at 
the right hand side of Eqs. (3.6) and (3.7) are: 



9' +ne^-v 



^ + ^^(^ + ^) 



0, (3.9) 



t-V^ 



*\' n^-n' 



n)^^{r-^ 



+ VXv - VV^ = 0. (3.10) 



The quantity Oi, appearing in Eq. (3.10) is the neutrino anisotropic stress. In fact, since the 
neutrinos are collisionless particles, their anisotropic stress (connected to the quadrupole moment 
of the neutrino phase space distribution as o^, = Ty^l'^ is not erased after neutrino decoupling 
occurring, approximately, when the plasma was dominated by radiation and for temperatures of 
the order of 1 MeV. The total anisotropic stress introduced in Eq. (2.35) will then be written, for 
convenience, as 

5,9^n;. + d^(P% = (p, + p,)vV, + (p^ + p^)vVb, (3.11) 
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where ub is the anisotropic stress of the magnetic fields. RecaUing the exphcit form of 11^ given in 
Eq. (2.23), it is easy to deduce that 



VVb = 



= 7 — 747 + 7^ ^B, 



a^p^ 



(3.12) 



where it has been used that V • J = and that AttJ = V x i?. In Eq. (3.12) the ratio bewteen the 
magnetic energy density and the photon energy density: 



^Bix) 






(3.13) 



Since p-y 



', Eq. (3.13) imphes that 0,b{x) is constant in time. Therefore, according to Eqs. 



(3.12) and (3.13), the contribution of the magnetohydrodynamical energy-momentum tensor to the 
scalar fluctuations reduces, effectively, to the magnetic energy (and pressure) densities (i.e. 6pB 
and (5/9b/3) and to the divergence of the Lorentz force. The collisionless nature of the neutrinos 
implies that not only the quadrupole moment (i.e. the anisotropic stress) of the neutrino phase 
space distribution contributes to the fluid evolution but also the higher moments should be taken 
into account. In particular, the following gauge- invariant equations have to be introduced 



cr„ 



:Fr 



4 

— ( 
15 



^/C.F.3, 



k 



-[^■^. 



u(e~i) 



+ '^)^u {£+!)]■ 



(3.14) 
(3.15) 



2£+l' 

Equation (3.15) applies for multipoles larger than the quadrupole. To specify initial conditions for 
neutrinos would then require, in principle, all the multipole moments of the (perturbed) neutrino 
phase space distribution. In practice, it will be sufficient to analyse the initial conditions of the lower 
moltipoles i.e. monopole (density constrast), dipole (peculiar velocity) and quadrupole (ansotropic 

stress). 

Unlike the CDM and neutrinos, the evolution of the baryons and of the photons is coupled by 
Thompson scattering: 



e'=v' 



n) ^^W 



e'^ + m^ = v^ 









V^C7 + aneXeaT(^b 
■(7^ + ^ 



^77' 



V • [ J X B] Apy 
Pha'^ 3 pb 



(3.16) 

^b), (3.17) 



where ut is the Thompson scattering cross section, n^ is the electron density and Xe is the ionization 
fraction. By taking the difference of Eqs. (3.16) and (3.17) the evolution of the combination (0^ — 9]^) 
can be obtained with the result that 



{6y — Oh)' + aUeXeCTT 



1 + 



4p^ 
3pb 



H^b - V2(C^ + XI,) 



V-[JxB] 



Pha 



(3.18) 
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For T <^ Toq the characteristic time scale for the synchronization of photons and baryons is set by 
(neUx)^^ which is smaU in comparison both with the expansion time, i.e. ar, and with the typical 
oscillation time, i.e. ar/k where k is the Fourier wavenumber. Consequently, any deviation of 
{6^ — 6\j) rapidly decays away. This aspect can be appreciated from Eq. (3.18) by taking the limit 
(Tt -^ oo. No matter how strong the forcing term is, (6^ — 6]^) -^ 0. Because of the tight coupling 
between photons and baryons, the baryon-photon velocity as a single physical quantity, i.e. 9^h 
whose evolution equation can be obtained by multiplying Eq. (3.16) by {4/3)p^ and by adding it to 
Eq. (3.17) (multiplied by pb)- The scattering terms then cancel and the evolution of ^^b) becomes: 



a/ "7b >/ __y^ 

(,7--—, '.b--^: 



Rh + l 



'7b + IT^T^^b - V 






V^i^ + C'y) 4VVB-V2f] 



B 



l + i?b + 4(l + iib) '^^'^^^ 

To derive the last term in Eq. (3.19), Eq. (3.12) has been used to express the Lorentz force in 
terms of the anisotropic stress and of the magnetic energy density in units of the photon energy 
density. In Eq. (3.19) the baryon to photon ratio has been also introduced as 

where z + 1 = l/a(r) is the redshift with the conventions adopted in Eq. (3.1). Since Eq. (3.19) 
involves only gauge-invariant quantities, it is simple to obtain each corresponding equation in a 
specific gauge. For instance, the conformally Newtonian gauge is often employed for numerical 
discussions of the tight-coupling approximation. In this gauge the perturbed metric components as 
simply 6sgoo = "^o^cf) and 6sgij = 2a^^(5ij. The gauge-invariant quantities appearing in Eq. (3.19) 
read, in the longitudinal gauge, 

C,^-« + ^, Ch^-i- + j, «,b^«,b, (3^21) 

*-*. R^-,.-^(,4) (3.22) 

where 5^ and 6^ are the density contrasts computed in the conformally Newtonian gauge. Equation 
(3.19) becomes, after easy algebra, 

,5; = 4^-' - ^e^b, <5(, = 3V''-e-yb, 



^^^^:^v^^'""i(m^"^'^+ 4(i+i?b) 



+ 17^^7b = -TTT^T - ^ <t^ + -TT^nr^^ (3-23) 



which is the standard form to the tight coupling equations in the conformally Newtonian gauge if 
cjb = ^B = 0. It is useful to mention that the total ( (see Eq. (2.33)) is given by 

Pt Pt Pt Pt 3(pt+Pt) 
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where p^-, P7, Pb and pc are, respectively, the unperturbed energy densities of neutrinos, photons, 
baryons and CDM particles. This observation is relevant for the actual solution of the Hamiltonian 
constraint (2.27). Finally, to solve the momentum constraint (2.19), we recall that the total peculiar 
velocity 9 is given, in terms of the peculiar velocities of each fluid, by the following expression 



(Pt + Pt)0 = {Pu + Pu)Ou + (Pf + P-y)0-y + (Pc + Pc)Oc + {Ph + Pb)^b- 

With the same notations of Eq. (3.24), p-y = p.y/3, py = Pu/'i while Ph = Pc = 0. 



(3.25) 



3.3 The adiabatic mode 

By going to Fourier space, Eqs. (3.6) and (3.7) can be easily solved, to lowest order in /cr < 1 and 
for r <C Teq. Denoting with 7?.* the curvature perturbation on super-Hubble scales and with $^, 
value of the Bardeen potential in the same epoch we will have that the Hamiltonian and momentum 
constraints imply 



^* = -2- 



5 + 2Ry 



R^{AaB-Ru^B) 



(3.26) 

(3.27) 

(3.28) 



15 + 4Ru 15 + 4Ry 

In Eq. (3.28) the fractions of neutrinos and photons present in the radiation-dominated plasma 
have been introduced according to the following conventions: 



R^ 



1 — Ry, 



Ri, 



1 + r' 



7 /4\4/3 

8 \llj 



(3.29) 



with Ni^ = 3, r = 0.681. The gauge-invariant density contrasts will then be determined to be 

^B 



C-, = C. = Cb = Cc = ^* - '-fR^. (3.30) 

By explicit integration of the corresponding equations, the expressions of the peculiar velocities are 



-A:V(7^* + ^,), ey = -k'^T 



, ^ ,„ N "B „ R'y 

4 K,j 



77 b 



-fcV 






Finally, the expression of the neutrino anisotropic stress turns out to be: 



Ov 



R^ feV 

Ry 3Ru 



n^ + -^, 



(3.31) 
(3.32) 

(3.33) 



From Eq. (3.28), ^* is determined from IZ^ and from the parameters of the magnetized plasma, 
therefore the ^^, appearing in the other formulae have been used just to avoid too cumbersome 
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expressions. In the limit 0,b -^ and ub ^ the solution reproduces the standard adiabatic 
mode. Using Eqs. (3.21) and (3.22) relating the gauge-invariant quantities with the conformally 
Newtonian fluctuations, we obtain, for instance 



■0* = </'* 



1 + Ir. 

5 



+ ^(4aB - Ru^b), TZ* 



^{15 + AR^) + ^{Aan 



RM 



iyi'Bj 



(3.34) 



If r^B -^ and 0"B -^ 0, Eq. reproduces the well known result with the difference of the Newtonian 
potentials being constant (to zeroth order in kr) and proportional to the fraction of neutrinos [46]. 

3.4 The non-adiabatic modes 

The discussion presented in the case of the adiabatic mode can be repeated, with the due differences, 
in the case of the isocurvature modes. Here only few examples will be given with particular attention 
to those features that will be useful for the forthcoming numerical aspects of the analysis. In the 
case of the neutrino isocurvature density mode, to lowest order in kr, the Bardeen potential is 
given by the following expression^ 



^ = ^, 



4 Ry Rp 



3{4Rp + 15) 



'3 3^ 

•S* - 75-0"B + Ti'B 
K,, 4 



(3.35) 



where the term S^, is the specific entropy of the photon-neutrino system. The term 5* is given, in 
the present approach, exactly by the mismatch between the density contrasts of the photons and 



of the neutrinos, i.e. S, 



1^7 



then: 



Cu 



-3(Czy — C7) = "S**- The density contrasts of each individual species are 

C-y = ^S,-^^B, Cc = Cb = -^^B. (3.36) 



^_1^5, + ^17b), Cy = ^S,-^^B, Cc = Cb = -%B. 



3 V 4 "7' " 3 i ' ' ' 4 

Again, integrating the evolution equations of the peculiar velocities the following results can be 
obtained: 

4 



-feV^. 



-fcV 



^s, + ^, + ^^ 



Rv 



O-B 



77b 



-feV 



2M/, + ?^S, 



i?7 + l^ 

— '- ilB + 0-B 



(3.37) 

(3.38) 



A feature of the neutrino isocurvature density mode is that the total density constrast, C, vanishes 
to zeroth order in /cr ^ 1. Consequently also the total curvature fluctuation, i.e. TZ, vanishes in 
for typical wavelengths larger than the Hubble radius. The neutrino anisotropic stress is given, in 
this case, by 



2^2 



(Tv 



Ry k^T 

Rp 2Ru 



** 



(3.39) 



As in the case of the magnetized adiabatic mode, the only free parameter of the solution is <S, ; the value of ^, 
is determined by <S, according to Eq. (3.35). For practical purposes, however, we will give some of the quantities as 
a function of both <S, and ^t,. 
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It is useful to make contact with the conformally Newtonian gauge where, according to Eqs. (3.21) 
and (3.22), the two Newtonian potentials are 



3{4R^ + 15) 



O^ 



R, 



3 3^ 

(Tb + -iiB 



1p* 



4:RjRi, 

'3(lR~fl5) 



*-5* 



R, 



3 3^ 

(JB + -i2B 



(3.40) 



Notice that 2^^, + (/)* = 0. This implies, in the conformally Newtonian gauge, that TZ = — (^*+(^^,/2) 
vanishes, as already pointed out by looking directly at the gauge-invariant fluctuations ( and 7^. 

Consider finally the situation where the non vanishing fluctuations in the specific entropy are 
given by 5c 7 and 5c jy and define, conventionally. 



-"07 



Sr 



o^. 



(3.41) 



This is the case of the CDM-radiation isocurvature mode. Then, from Eq. (2.31), by making the 
sum over the fluid components explicit we do obtain 



. _ Pcipy + P'u) e 

9Hpt 



(3.42) 



where we used that S^c = Syc = —5*. Now defining Pr = P7 + Pu-, Eq- (3.42) can also be written 



as 



SPn&d = Pc Cg 5* 



Pr 



4a. 



eq 



(3.43) 



' 3 3/Jc + 4pr 3(3a + 4aeq) ' 

where the second equality in the definition of the total sound speed Cg follows from the conventions 
illustrated in connection with Eq. (3.1). Therefore, to lowest order in /cr < 1 and for r < Teq we 
have that 



7^~C 



45, + SQbR^ 



16 



'■cq 



* 



15 + 4Ry 
;(15 + 2Ry 



5* + —^^bR"/ 



-"■eq 



(3.44) 
(3.45) 



where we used that fi^ <^ f^c (so that acq — ilf/ilc). In the limit f^B -^ 0, Eq. (3.44) implies 
that TZ ~ — 5*(a/acq)/4 which reproduces the standard result obtainable for the CDM-radiation 
isocurvature mode in the absence of magnetic fields. In the same limit of vanishing magnetic con- 
tribution (and with vanishing fractional contribution of neutrinos, i.e. R^ ^ 0) ^ ^ 5^,(a/acq)/8 
which is also well known from the analysis of the two fluid system in the absence of neutrinos. In 
this case the relevant density contrasts are given by 



C7 — Qf 



Qt 



-Ry, 



Cc 



"»_>=): 



Qt 



'-Ry, 



(3.46) 



as it follows from the Hamiltonian constraint and from the definition of ( in terms of the density 
contrasts of each fluid. The peculiar velocities can then be directly obtained to lowest order. In 
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particular it is interesting to comment about the neutrino sector where the solution for the lowest 
multipoles can be written as 



(Jv 



15 Qc 
y 15 + 2R^ 



Sif + —0,BR-y 



fc^r^ + kr 






2^3 



fcV 



n. 



R 



•^^^-15 



R^ 
Rv 



-CJB 



6 15 + 2i?^ 



Sif + —Q.'qR^ 



kr. 



(3.47) 



Notice that to set correctly initial conditions in the magnetized case the octupole is important. 
However, in the limit of vanishing magnetic inhomogeneities, the octupole can be simply set to 
zero as all the remaining terms in the hierarchy. Finally, in the case of the CDM-radiation mode, 
using Eqs. (3.21) and (3.22), the conformally Newtonian potentials are: 



^ = 8(15 + 2R^) y^* + 4^^^^ 



-"■eq 



15 - AR^ 
8(15 + 2R^ 



s^ + —vi^qR^ 



^eq 



(3.48) 



The two longitudinal potentials V and (j) are only equal in the limit Ri^ -^ 0. In the two-fluid 
treatment of the problem this would actually be the case since no neutrinos are allowed. 

4 Magnetized curvature fluctuations after equality 



The regular solutions describing both adiabatic and entropy modes will now be followed through the 
transition between the radiation-dominated stage of expansion and the matter-dominated epoch. 
In fact, as already mentioned, the decoupling time takes place after Teq (i.e. r^cc 
h^Tirn = 0.134). 



2.36 Tcq for 



4.1 The adiabatic mode 

In the absence of magnetic fields and in the absence of non-adiabatic pressure density fluctuations, 
Eqs. (2.27) and (2.29) imply that the total curvature perturbation TZ, i.e. the curvature fluctuation 
on comoving orthogonal hypersurfaces, is conserved for wavelengths much larger than the Hubble 
radius, i.e. for kr -^ 1. Conversely, the Bardeen potential ^ diminishes, in the transition from 
radiation to matter of a factor 9/10 (assuming vanishing neutrino anisotropic stress). In the absence 
of magnetic fields, the effect of the neutrino anisotropic stress is often ignored for order of magnitude 
estimates. However, if the plasma is magnetized, this approximation is no longer justified as 
clearly demonstrated by working out the solutions presented in section 3. In fact, the anisotropic 
stress of the neutrinos is balanced by the magnetic anisotropic stress and the subsequent evolution 
of the system will tell what happens once the initial conditions are accurately set. Adiabatic 



19 



(magnetized) initial conditions are imposed according to Eqs. (3.26), (3.27) and (3.28). The 
initial conditions for the density contrasts can be found in Eq. (3.30). The velocity fields and the 
neutrino anisotropic stress are then the ones derived, respectively, in Eqs. (3.31)-(3.32) and (3.33). 
By setting JIb = ctb = in the initial conditions the standard adiabatic mode (with canonical 
dependence upon the fractional contribution of massless neutrinos) is reproduced. 

In the realistic situation the adiabatic mode will have a nearly scale-invariant spectrum (with 
red tilt) while ^b will have a nearly scale- invariant spectrum (with blue tilt). For simplicity the 
relevant results concerning the power spectra are collected in the appendix here we just recall that 
the spectra of 7^* and ile and ctb can be written, according with the notations of the appendix, as 

/ k \ "'■"^ 
Vn(.k)=An[-j , (4.1) 

and as 

_ ^^^ ^. ^7=T2 /A;\2^ „,,, (3-2e)(188-4e2_66e)_ ,,, ,, , 

where kp = 0.05 Mpc~^ is a conventional pivot scale of the scalar adiabatic modes and ki = \ Mpc~^ 
is the pivot scale usually employed to parametrize the magnetic power spectra. The quantity e is 
the spectral tilt of the magnetic energy spectrum i.e. < e < 1. Through the years it emerged 
that a nearly scale invariant magnetic energy spectrum is probably the most interesting to discuss 
in CMB studies since it is not strongly affected by diffusive processes (see, for instance, [10, 12]). 
The other numerical quantities appearing in Eq. (4.2) are defined in Eqs. (A. 15) and (A. 16) of 
the appendix. If the magnetic fields are absent and n^ = 1, the Sachs- Wolfe plateau (see section 5) 
corresponds to 

^^Q = 4f, ^^ = 2.65x10-^ (4.3) 

ivr 25 

The first interesting exercise (even if not realistic as we shall see later on) is to assume that the 
magnetic field is exactly comparable in magnitude with the curvature perturbations, i.e. IZ^, = VL-q- 
In particular without loss of generality, we will assume rir = 1 and e = 0.01. In this situation, 
as dictated by Eq. (3.33) total anisotropic stress S(t) = a^ + Rj/RuO'b is initially vanishing for 
typical wavelengths larger than the Hubble radius. Figure 1 (plot at the left hand side) shows that 
the total anisotropic stress remains vanishing as long as the mode is outside the horizon and it 
oscillates as soon as the mode reenters. In Fig. 1 the time around the radiation-matter transition 
is measured in Mpc while the comoving wave-numbers are measured in Mpc~ . In the right plot of 
Fig. 1 we use the same initial conditions but with a value of f^B that is one tenth of the amplitude 
of the initial curvature perturbations 7^*: the amplitude of the oscillations of the total anisotropic 
stress diminishes. The evolution of the curvature perturbation TZ presents some interesting aspects 
that allows a comparison with the analytical results. In Fig. 2 (left plot), the evolution of TZ is 
reported for different wave-numbers and for the same range of parameters chosen in the case of 
Fig. l(left plot), i.e. f^B = 7^*. When the given wavelength is larger than the Hubble radius during 
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Figure 1: The time evolution of the total anisotropic stress for two different amplitudes of the 
magnetic field intensity. The comoving wave-numbers are measured in Mpc^ . In the plot at the 
left J7b = 7^*. In the plot at the right Ob = 0.17^*. 

the radiation epoch(i.e. r <^ Tcq), C' — TZ' — 0. This conclusion can be derived, for instance, 
from Eq. (2.29)) and from the Hamiltonian constraint (2.27). Equation (2.29) tells that when 
Cg ~ 1/3 (and when 5pnad = 0, as in the case of the adiabatic mode), C, (and therefore TZ, by virtue 
of Eq. (2.27) ) are conserved for kr <^ 1. As soon as the given wavenumber reenters the Hubble 
radius the conservation of TZ no longer holds and this can be appreciated from Fig. 2 (left plot). 
The interesting situation occurs when the given wavelength is larger than the Hubble radius after 
equality (but before decoupling). These are the modes determining the large-scale contribution to 
the (scalar) CMB anisotropics and, in particular to the SW plateau. From Fig. 2 it is clear that the 
magnetic field drives a growth of TZ when the given mode is still outside the horizon. In the right 
plot of Fig. 2 the evolution of TZ is compared with the case where the magnetic field is absent (the 
numerical values of the other initial conditions are the same) . The horizontal line with amplitude 
1 is the value of TZ in units of TZ^ when the magnetic contribution is absent. The middle curve 
represents the contribution of the magnetic field when Qb = o"b = 7^* • The dashed line in the right 
plot of Fig. 2 is the analytical result obtained by solving Eqs. (2.29) and (2.27) for kr <C 1. The 
analytical result reproduces very well with the numerical estimate. Consider, indeed, Eq. (2.29) 
recalling that within the notations introduced so far 



''" 3(3a + 4aeq) ' 
Then Eq. (2.29) becomes, in the limit fcr ^ 1 

da 



w{a) 



-''eq 



3(a + Oeq) 



R-yil'QCLl 



eq 



(3a + 4aeq)2 ■ 



(4.4) 



(4.5) 



The initial conditions deep in the radiation epoch (i.e. a/ceq — > 0) are such that the Hamiltonian 
constraint must be satisfied, i.e. C* — T^* + ©(Aj^t^). Consequently, the solution of Eq. (4.5) and 
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Figure 2: The time evolution of curvature perturbations for H.^ = IZ^. The dashed hne in the right 
plot is the analytical result valid for wavelengths larger than the Hubble radius. In both plots TZ is 
rescaled by its initial value, i.e. 7^*. 



(2.27) will be 



7^ ~ c{a) = n^. 



4(3a + 4aeq) ' 



In the limit a/oeq > 1 the above expression leads to 



7^~7^, 



R^ 



2_2\ 



J^B + 0(fc r 



(4.6) 



(4.7) 



For J7b — Tl* we get that TZ/IZ^ ~ (1 — Rj/4:) ~ 0.85 which is the asymptotic value reported, 
with the dashed line, in Fig. 2. The same result can be readily obtained from Eq. (2.36) since 
the neutrino anisotropic stress counterbalances the magnetic stress with the net result, as verified 
numerically in Fig. 1, that the total stress can indeed be neglected for kr <^ 1. Once TZ and (^ have 
been computed, ^ can be derived from Eq. (2.35) whose explicit form can be written as 



d^ 5a + 6 ^ 

+ 1^-, —'i' 



3f^B 



7^*■ 



3a + 4 



da 2a(a + l) 8(a + 1) 2a(a + 1) 

where the variable a = a/ogq has been introduced. Recalling now that 



din /(a) 5a + 6 



da 



fia) 



a 



Eq. (4.8) can be written as 



d 



— [/(a)^(a)] 
da 



2a(a + 1) 
^_^^a'(3a + 4) 



+ 



VaTT' 



snsa^ 



'2(a + 1)3/2 8(a + 1)3/2' 
Direct integration of Eq. (4.10) between and a implies 



(4.8) 



(4.9) 



(4.10) 



fia)^ia) - /(O)^(O) 



l6{VaTT - 1) + 9a3 + 2a2 - 8a 



+ 



-7^* ^ 

15Va + l 

3R^nB 32(1 - \/^^TI) + 2a3 - 4a2 



16a 



40 



Va + 1 



(4.11) 
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Now, since /(O) goes to zero quite fast for a ^ 0, the term /(O)^(O) -^ provided ^(0) does not 
diverge in a ^ 0. Therefore Eq. (4.11) imphes 

16(Van - 1) + 9a^ + 2a^ - 8a 



^(a) = -7^*- 



15a3 



SR^ftB 32(1 - V^TT) + 2a[8 + a(a - 2)] 

+ ^n q • (4.i2j 

In the limit a > Ocq (i-e. a > 1), Eq. (4.12) implies that ^(q) goes to a constant (in time) and 
that, in particular, 

liui^{a)^-ln, + ^nB. (4.13) 

Q>i 5 20 

The reported results are gauge- independent, this means that if we would perform the numerical 
integration in a specific gauge (hke, for instance, the conformally Newtonian gauge) exactly the 
same results for TZ and C will be obtained. This exercise has been done. We integrated explicitly 
the evolution equations in the longitudinal gauge and obtained exactly the same results derived by 
numerical integration of the gauge-invariant system of equations. 

As an example, we can illustrate the equivalence of these two approaches by presenting the 
numerical evaluation of the Sachs- Wolfe contribution in the sudden decoupling limit where the 
visibility function is approximated by a Dirac delta functions 6{t — Tdcc) and the optical depth is 
simply its integral, i.e. a step function 9{t — Tdec)- In this (naive but physically useful) limit the 
brightness perturbation (related to the Stokes parameter I) can be written as 

/•To 

Ai(fc,ro)= / e^'=^(^-^")[V' + <^']dT + e^'=^(^<^--^«)[Aio + + MF],,^^, (4.14) 

where fi = k -h is the angle between the direction defined by the momentum of the photon and the 
direction defined by the Fourier mode k. The term Ajq is exactly 5^/4 while 

The first equality in Eq. (4.15) is just a redefinition; the second equality holds in the tight coupling 
approximation. In Fig. 3 the result of the numerical integration is illustrated for the case where 
Ob = TZ^:. The final time in the linear scale appearing in both plots of Fig. 3 is Tdcc (recall that 
Tdec = 282.81 Mpc and r^q = 119.07 for the typical values of the parameters discussed in Eqs. (3.2), 
(3.4) and (3.5)). In the left plot of Fig. 3 we illustrate the combination A(r) = 6^/4 + (j) (recall 
Eqs. (3.21) and (3.22) for the definitions of the conformally Newtonian fluctuations). In the plot 
at the right of Fig. 3 we illustrate V{t). 

The results of Fig. 3 should be compared with the ones of Fig. 4 where the same quantities are 

illustrated but for JIb = 0.001 IZ^. The comparison of the two plots shows that the introduction of 

^This is not only the case for the magnetized adiabatic mode but also for the majority of the isocurvature mode 
with the exception of the neutrino isocurvature velocity mode. In fact, in the case of the adiabatic mode Vl/(0) is 
constant but /(O) -^ (sec Eq. (4.9)). In the case of the isocurvature modes ^(0) -^ 0. 
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Figure 3: The evolution of the Sachs- Wolfe terms for the case Qb = Tt* and for three typical 
wave-numbers in units of Mpc~ . 
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Figure 4: The same quantities illustrated in Fig. 3 but now for Ob = 0.0017^* 
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Figure 5: The evolution of the total anisotropic stress in the case of the CDM-radiation (magnetized) 
non-adiabatic mode. In the plot at the left Ob =5*. In the plot at the right f^B = 0.15*. 

magnetized initial conditions affects the amplitude of the oscillations and, less strongly, the phases. 
For Qb = the results coincide with the ones of Fig. 4. This means that in order not to affect 
the dominant adiabatic mode, JIb should not exceed 10~'^7^*. We will get back to this point in 
section 5. It is appropriate to note here that the oscillations appearing in Figs. 3 and 4 are the 
Sakharov oscillations that will determine the Doppler peak structure. The requirement that the 
dominant adiabatic mode is not affected by the presence of the magnetic field intensity implies then 
a standard Doppler peak structure. 

4.2 The isocurvature mode 

Equations (4.5) and (4.6) suggest that the magnetized contribution can indeed be seen as an 
isocurvature correction to an observationally dominant adiabatic contribution. This aspect can be 
appreciated by looking at the second term at the right hand side of Eq. (4.6): notice that this 
term goes to zero as a/oeq <C 1. It is therefore natural to consider an example of a magnetized 
isocurvature mode, like, for instance, the CDM-radiation isocurvature mode. In this case the proper 
(magnetized) initial conditions for the numerical integration can be found in Eqs. (3.43)-(3.47). 
Notice, again, that by setting ub ^ and I^b ^ we do recover the initial conditions for the 
standard CDM-radiation isocurvature mode with non-vanishing fractional contribution of massless 
neutrinos. In the case where J^b = 5* (where 5* parametrizes the primordial entropy fluctuation 
in the CDM-radiation system) the anisotropic stress evolves as in Fig. 5 (left plot). Conversely, if 
Ob = 0.15*, the numerical results are reported, always in terms of the total anisotropic stress, in 
the right plot of Fig. 5. The two plots present the same features already discussed in the case of the 
adiabatic mode (see Fig. 1). In Fig. 6, for comparison, the behaviour of S(rdec) is illustrated as a 
function of the wave- number in the case of adiabatic initial conditions (left plot), and in the case 
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Figure 6: The total anisotropic stress is illustrated at a fixed time (coinciding with Tdec) but as a 
function of the wave-number measured in Mpc^ . In the left plot the adiabatic case is illustrated. 
The right plot illustrates instead the CDM-radiation mode. In both cases different strengths of the 
magnetic field intensity have been reported. 



of isocurvature initial conditions (right plot). The different values of the magnetic field strengths 
are referred to the amplitude of the initial (adiabatic or non-adiabatic) mode. 

The crucial difference between the (magnetized) adiabatic mode and the (magnetized) CDM- 
radiation mode can be illustrated by looking at the evolution of 7^ and ^. In contrast with the 
adiabatic result, however, both ^ and Tt vanish for a/a^c^ — > 0. This feature is quantitatively 
illustrated in the plots of Fig. 7. In the right plot the evolution of ^ is reported in the asymptotic 
region where the analytical estimates are easily compared with the numerical results. Both for IZ 
and ^ for different sets of magnetized isocurvature initial conditions the analytical results reproduce 
correctly the numerical ones. The analytical derivation resembles, in this case, the one already 
performed in the adiabatic case with the difference that, now, the non-adiabatic pressure density 
variation is given by 

[Pt + Pt) (4aeq + Say 

By means of this expression, Eqs. (2.29) and (2.27) imply that, for a/ccq ^ 1 



3 4 



5 20 



(4.17) 



The results obtained so far in the case of the adiabatic and in the case of the isocurvature mode 
can be summarized by the following two matrices, i.e. 

'7^*■ 




'Mnn Mns Mub' 
Mss MsB 
Mbb 



(4.18) 



26 



R(t) 



-0.2 
-0.3 
-0.4 
-0.5 





X 






"^ 






\ 




/ 


k=o.o: 






















/ 








k=0. 


DOl\\^ 




k=0.0001 



2 3 4 5 

Log(T/Mpc) 



0.3 

0.25 

0.2 

i(r) 0.15 

0.1 

0.05 



















7 i 

- ■' 1 










\ k= 0.001 

1 






/ 






1 
k=0.C001 




/ 












/ 






k=0.0 


L 


y 








<^ 





2 3 4 

Log(i/Mpc) 



Figure 7: The evolution of the magnetized CDM-radiation isocurvature mode is ihustrated in terms 
of TZ and ^ for JIb = 5* . 

In the case of a mixture of (magnetized) adiabatic and CDM-radiation modes, we find, for a > Oeq 

1 . . R-. 



M 



Tin 



1, 



M 



ns 



Mtzb - ^, 



M 



ss 



1, 



M 



SB 



0, 



(4.19) 



3 ' '"^4 

and A^BB -^ 1- The matrix (4.18) gives the values of the various fluctuations after equality in 
terms of the initial conditions set for r <C Tgq. In similar terms, for ^ we will have 





V 




where, now, 



A^*-^ -^ TTT 



9 

fo' 



Mq,S -^ -: 

5 



M 



*B 



oRy 



(4.20) 



(4.21) 



while the other matrix elements are the same as in the case of Eq. (4.18). In Fig. 8 the curvature 
perturbations are illustrated at T^ec but as a function of k, i.e. the comoving wave-number. For 
wavelengths larger than the Hubble radius, the flat plateau appearing in both plots of Fig. 8 is 
accurately predicted by Eqs. (4.18) and (4.19). 

The techniques presented here may have useful implications. For instance, in [47] it has been 
correctly argued that magnetic fields may be a source of non-Gaussian signatures in CMB maps. 
It can be easily argued that if the magnetic fields are Gaussianly distributed, the magnetic energy 
density will be non-Gaussian (being the square of a Gaussian variable). This argument is usually 
made by looking at the magnetic energy-momentum tensor. Here it has been shown that , in the 
long wavelength limit, the magnetic contribution to curvature perturbations is exactly proportional 
to the magnetic energy density with a numerical coefficient that is dynamically determined by the 
features of the radiation-matter transition. The results of [47] on higher order correlators must then 
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Figure 8: The gauge- invariant curvature perturbations TZ are reported at fixed time (coinciding 
with Tdec) and as a function of the comoving wavenumber in units of Mpc~^. In the left plot the 
adiabatic case is illustrated. In the right plot non-adiabatic initial conditions are imposed. 

be complemented by a consistent framework where adiabatic and non-adiabatic initial conditions 
are properly set. 

It would be interesting to go to higher order in the tight-coupling expansion and to adapt the 
semi-analytical techniques already developed in the standard case [49, 50]. This will allow a semi- 
aanalytical treatment of the implications of magnetized initial conditions on polarization power 
spectra. 

5 The Sachs- Wolfe plateau for magnetized initial conditions 

Magnetized curvature and entropy perturbations, initially present for r <C Teq, have been followed 
through the radiation matter transition. This analysis allows the calculation of the observables 
that determine the large-scale temperature fluctuations when the relevant modes are outside the 
horizon after equality but before decoupling. While this calculation is much simpler in the absence 
of magnetic fields, it has been shown that the inclusion of magnetic fields well defined since the 
regular modes present deep in the radiation epoch propagate regularly through Teq. It is essential 
to treat carefully not only the magnetic anisotropic stress but also the anisotropic stress of the 
fluid (mainly induced by collisionless neutrinos). Therefore, a posteriori, it is also clear why, in 
the absence of magnetic fields, the analytical calculation is much simpler: in the standard case the 
neutrino anisotropic stress can be neglected and its inclusion only produces a weak dependence on 
the neutrino fractional contribution that enters, for instance, in the relation between the curvature 
perturbations 7^ and the Bardeen potential. It is then appropriate to proceed further and to eval- 
uate the Sachs- Wolfe (SW) plateau induced by magnetized curvature perturbations. This analysis 
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may serve as a starting point for future numerical investigations aimed at the inclusion of regular 
magnetized modes in the discussion of parameter extraction from present and forthcoming CMB 
data. The temperature fluctuations induced by the scalar modes, i.e. the brightness perturbations 
of the intensity of the radiation field, can be written as (see, for instance, [28]) 

Ai(f , n, r) = — ^ I d'ke-^'^-^Ai{k, h, r), (5.1) 

where n represents the direction of the photon momentum and, for convenience we also introduce 
^ = h ■ k, i.e. the projection of the photon momentum on the direction of a given Fourier mode. 
Notice that the subscript I appearing in Eq. (5.1) refers to the first Stokes parameter measuring 
the intensity of the radiation field. The Fourier amplitude appearing in Eq. (5.1) can then be 
expanded in a series of Legendre polynomials as ^ 

Ai(^, n, r) = ^(-i)'(2^ + l)Au{k, T)Pe{fi). (5.2) 

e 

Assuming then that the receiver is located at x = 0, the temperature anisotropics can be expanded 
in spherical harmonics as 

Ai(n, r) = ^a^„Y£„(n). (5.3) 

£ m 

Hence the coefficients a£ m can be determined as 

aem = -7^;^^{-i)^ I d^kYi^Au{k,T), (5.4) 

where the theorem of addition of spherical harmonics [51, 52] has been used. 

Since the magnetic field is fully inhomogeneous the isotropy of space-time is preserved and, 
therefore, the angular power spectrum will not depend on the azimuthal angle. In particular the 
magnetic field is described in terms of a two point correlation function (see Eqs. (A. 5) and (A. 6)) 
that, therefore, does not imply any specific preferred direction. This means that the ensemble 
average of the ai m can be expressed as 

{aema}'m') = Ci6ii'6mm'- (5.5) 

Using again the addition theorem for spherical harmonics [51, 52] the angular power spectrum will 
then be expressed directly in terms of the Cg, i.e. 

C{6) = (Ai(ni,ro)Ai(n2,ro)) = ^ ^(2£ + l)QP,(m • ns). (5.6) 



*We report the following formula in order to set accurately the notations. In fact different authors choose to 
perform the deceomposition in terms of Legendre polynomials following different conventions. For instance, some 
authors include the factors (— J)* (or even (2i? + 1)) in the definition of A^. Of course these difi'erences are immaterial 
once the calculation is performed consistently in the same framework. 
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We have now to compute Ci in terms of the curvature fluctuations. In order to do so we solve the 
heat transfer equation in the sudden decoupling limit. The large-scale contribution ^ to Ai(fe, r) at 
the present time tq can then be written as 



Ai(/c,ro) = e 



-ifc^(ro-rdec) 






(5.7) 



where Tdcc ^ ''"o and its presence can be neglected. The terms at the right-hand side of Eq. (5.7) 
represent the contribution of a general perturbation to the temperature anisotropy. 

If S-t: = r^B = the only contribution to the angular power spectrum comes from the usual 
adiabatic mode. The dependence upon the primordial perturbation is also canonical. For instance, 
recalling that, in the longitudinal gauge, TZ^ = — (5/3)'(/'m we have that the contribution is V'nf/3 
where V'j^ is the constant adiabatic mode during matter and in the conformally Newtonian gauge. 
Similarly, if the the only non-vanishing contribution is of isocurvature type we will have that 
■0m° = 5*/5 reproducing, therefore, the isocurvature contribution to the SW plateau that goes, in 
the conformally Newtonian gauge, as 2^J^°. We can then expand the plane wave appearing in Eq. 
(5.7) in Legendre polynomials, with the result, that Aj^ is determined to be: 



n^{k) 2^ ,-., 3ily^ ,^, 



^u{k,To) = 

where je{kTo) are the spherical Bessel functions i.e. [51, 52] 

je{kTo 



JeikTo), 



IT 



2kTQ 



Je+i/2{kTo). 



(5.8) 



(5.9) 



Now the calculation goes as follows. Equation (5.8) has to be inserted into Eq. (5.4). Then, the 
obtained result has to be inserted into Eq. (5.5) and the Cg can be determined from the explicit 
expression of the ensemble average of the aim- To compute explicitly {a^^apm') we have to know 
the ensemble averages of the various random fields appearing in Eq. (5.8). In particular, inserting 
Eq. (5.8) into Eq. (5.5) and trying to evaluate the expectation value (a^mO^'m) the following 
expectation value has to be computed 



{Au{k,To)Aii:{p,To)*) = 
9 

+ Y^{^Bik,To)nBip,To 



Je{kTo)ji>{pTo 
4 



1 - 4 - 

— (7^*(A;,ro)7^*(p,ro)*) + — (5*(fc,ro)5*(p,ro) 
25 25 



25 



(7^*(/c,ro)5*(p,ro)*) 



-^{^B{k,To)n^{p,To)*) + —{S^{k,To)QB{P,Toy 



(5.10) 



The stochastic averages appearing in Eq. (5.10) can be obtained by using the expressions appearing 
in Eqs. (A.2)-(A.4), (A. 13) and (A. 22). Recalling then Eq. (5.5) the expression for the Q is 



Ce = 4TT dlnkg{k)j]{kTo), 



(5.11) 



^We neglect here the integrated Sachs- Wolfe contribution as well as, for the moment, the Doppler contribution 
that is relevant for smaller angular scales 
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0{k) = ^^^^ + ^'Ps{k,ro) + -^R'^VB{k,To) - ±^Vnik,To)Vsik,To)cosyr 



where 



-■^^■yy'PBik, To)'Pnik, To) cos Jbr + W^Rjy'Psik, To)'PBik, To) COS7sb. (5.12) 

Concerning Eq. (5.12) the following comments are in order: 

• if the adiabatic, non-adiabatic and magnetic contributions are all uncorrelated (i.e. ^rs = 
7br = 7sb = "^72), the Sachs- Wolfe contribution will be enhanced if compared with the case 
where only the adiabatic mode is present; 

• if the adiabatic and the non-adiabatic modes are positively correlated (i.e. cos jrs > 0) but 
uncorrelated with the magnetic contribution (i.e. ji,r = "yts = "^72) the correlation reduces 
the contribution to the SW plateau; 

• if the adiabatic and non-adiabatic contribution are both uncorrelated (i.e. 'jrs = 7i"/2) but 
the adiabatic and magnetic components are positively correlated (i.e. cos7j.s > 0) again the 
amplitude of the SW plateau is reduced; on the contrary, if the non-adiabatic and magnetic 
components are positively correlated (i.e. cos 7^^ > 0) the amplitude of the SW plateau is 
enhanced; 

• the previous cases are of course reversed when the various modes instead of being correlated 
are anti-correlated (i.e. cos 7 < where 7 stands for one or more corrrelation angles). 

It is then clear that the presence of a regular magnetized mode during the radiation dominated 
epoch induces, after Teq super-horizon fluctuations whose effect on the SW plateau depend on the 
mutual correlation of the various components. This situation is, in some way, similar to the case 
where abiabatic and isocurvature modes are correlated (or anti-correlated). The present situation 
is, however, phenomenologically more rich due to the presence of the magnetic field. 

By performing the integral over k and by recalling Eq. (5.9), Eq. (5.12) allows the explicit 
evaluation of the Ce once the form of the power spectrum is specified. Each of the power spectra 
appearing in Eq. (5.12) will have a different spectral slope. So the general notation adopted here 
will be the one summarized in Eqs. (A. 3) and (A. 4) (for the autocorrelation of the adiabatic and 
isocurvature modes and their correlation), in Eqs. (A. 13) and (A. 22) (for the autocorrelation of the 
magnetic component and for the cross correlations between the magnetic sector and the adiabatic 
or non-adiabatic contributions). As discussed in the appendix, the power spectra of the various 
contributions are all referred to the same (pivot) wave-number kp whose specific value may range, 
according to present conventional choices between 0.05 Mpc~ and 0.002 Mpc~ . 
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With these necessary specifications the integration over k in Eq. (5.12) gives 
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Ci 



An 



9 



?2o2 



Zi{nr,ti + —AsZi{ns,t) + -r7^Rl^BL^2{(^J) - — VATzAsZi{nr s, i) COS Jr 



3 _ g _ 

— \/Ar,Rt ^B L2:3{nr,e,e) COS Jbr + —VAs Rj^B L ^s^Us, e, £) COS -fbs 



(5.13) 



where the functions Zi{m,i), Z2{e,i) and Z3{m,e,i) are defined, for a generic spectral index n 



z,{n,e) 
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fcn^^^ 



r(2-2e)r(^ + e) 

r^(l-e)r{e + 2-e) 



Vr ^ n + l 

— 2^2^~ 
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.F(e) 



fcL 



^/fcox^ni-^-tin^+i + f-i 



r2fZ_£_I»;]rf9_i_/'_e_n 
^'4 2 4l^l4^^ 2 4 



(5.14) 
(5.15) 

(5.16) 



The generic spectral index n appearing in Eqs. (5.14), (5.15) and (5.16) may correspond either to 
Ur, or to Ug or even to Urs = (f^r + ns)/2. Furthermore, the function .F(e) appearing in Eqs. (5.14) 
and (5.15) is given in Eq. (A. 15) of the appendix. 

As explained after Eq. (5.12), the phenomenological situation implied by the magnetized adia- 
batic and isocurvature modes is rather rich since we essentially have three spectral indices (i.e. n^, 
Us and e), three amplitudes (i.e. A-jz, As and ^bl), and three correlation angles^^. 

In order to discuss this situation in fully consistent terms it is necessary to perform a numerical 
analysis of parameters extraction similar to the one reported in [8] where, however, the magnetized 
contribution is taken into account both at the level of the initial conditions and at the level of the 
normalization. 

To get an estimate of the possible bounds obtainable in the present framework on the intensity 
of the magnetic field the minimal situation will be considered. Suppose, indeed, that the SW 
plateau is just produced by a magnetized adiabatic mode. In this case Eq. (5.13) reduces to 



Ci 



j[j^ 9 2 3 

—-Zi{nr,£) + -— R Sl-Q ^Z2{e, £) - — V^i?7 ^bl ^3("-r,e,^) C0S7 
25 lUU ' 25 



(5.17) 



^°The integrals in Eq. (5.12) converge provided —3 < n^ < 3, —3 < n^ < 3 and ~3 < n^s < 3. 
In principle the spectral slope of the cross- correlations may not be simply related to the spectral indices of the 
autocorrelation. In this case the formulae given in the present section are slightly modified. The modifications can 
be easily derived from the results of the appendix by assuming that, for instance, rir s is not given, as assumed here, 

by {ur +ns)/2. 
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Figure 9: The bounds on the magnetic field intensity derived from Eq. (5.17) assuming maximal 
anticorrelation between the magnetic and adiabatic components. 

where the correlation angle has been simply denoted by 7, since it is, in this case, the only correlation 
angle. Since the contribution to the SW plateau must be predominantly adiabatic we do know that 
the scalar spectral index must be nearly scale invariant. In the case where n^ = 1 and in the absence 
of magnetic field, Eq. (5.17) implies that i{i + 1)Ci/2tt = Ati/25. The WMAP determination will 
then impose An — 2.65 x 10^^. The recent release of WMAP data [3] seems to prefer slightly red 
spectra with preferred value (taking WMAP alone) rir = 0.951. The spectrum of the magnetized 
contribution is controlled now by e. If < e < 1 we are in the interesting situation where also 
the magnetic field spectrum is nearly scale-invariant (though with a blue tilt since e is directly 
related with the tilt of the magnetic energy spectrum Vn, in the notations of the appendix). The 
range < e < 1 is phenomenologically interesting not only because, in this range, the magnetic 
energy spectrum is nearly scale-invariant but also because the diffusive effects (i.e. finite thermal 
and magnetic diffusivity scales) do not affect the properties of the large-scale magnetic field (see, 
for instance, [9]). 

If the contribution to the SW plateau is predominantly adiabatic, then, the limit on the magnetic 
field intensity for a given range of rir and e are illustrated in Fig. 9 assuming the validity of Eq. 
(5.17) and enforcing the WMAP (i.e. experimental) determination of the SW plateau. In the left 
plot of Fig. 9 the maximal magnetic field intensity is illustrated as a function of the scalar spectral 
index rir for different values of e. For the selected values of the parameters the magnetic field has 
to be always below each of the curves. As discussed above (and also in the appendix), Bl is the 
magnetic field regularized over a typical comoving wave-number k^ that we take to be, according 
to the conventions, 1 Mpc~ . The pivot scale has been fixed at kp = 0.05 Mpc~ . To obtain the 
estimates illustrated in Fig. 9 we assumed, in Eq. (5.17) the less favourable situation, i.e. the 
one where 7 = vr. In this case (maximal antcorrelation) the two autocorrelations and the cross- 
correlation term all add-up to the plateau and there is no possibility of (artificial) reduction of the 
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amplitude. From the left plot of Fig. 9 it is clear that the present value of the magnetic field is 
bounded to be Bl < 10^^ — 10^^ G. Similar conclusion can be drawn from the plot at the right hand 
side of Fig. 9 where the allowed region is illustrated for a specific value rir (preferred by WMAP 
data alone) as a function of e. It is clear that by allowing a positive correlation (i.e. cos 7 > 0) 
instead of a negative one (as previously assumed) the bounds may be relaxed. It could be argued 
that the limit on magnetic fields arising from this analysis is competitive with the ones obtainable, 
for instance, from Faraday rotation of the E-modes into B-modes. The main conclusion is however, 
that it is now possible to include the magnetic fields in a systematic strategy of parameter extraction 
from available CMB data possibly combined with large scale structure data. 

6 Discussion 

There are no compelling reasons why the scalar modes induced by large-scale magnetic fields should 
not be considered as a possible source of CMB anisotropies. In this paper a new framework for this 
analysis has been proposed. The main findings can be summarized as follows: 

• the presence of large-scale magnetic fields (with nearly scale-invariant magnetic energy spec- 
trum) affects both the adiabatic and the entropy modes; 

• a general gauge-invariant system for the analysis of the initial conditions has been proposed: in 
such a system the gravitational sector is described in terms of the curvature perturbations on 
comoving orthogonal hypersurfaces, the density constrast on uniform density hyeper surf aces 
and the Bardeen potential; 

• this gauge-invariant system has been analytically solved deep in the radiation dominated 
epoch (well before equality) and then numerically followed after equality (but before recom- 
binantion) in the tight coupling approximation; 

• the numerical estimates (corroborated by the analytical results) allow to compute the pre- 
cise value of the magnetized adiabatic and non-adiabatic curvature perturbations when the 
relevant modes are outside the horizon before decoupling; 

• as an example the effect of magnetic fields on the SW plateau has been discussed. 

The magnetized curvature and entropy perturbations propagate then smoothly through equality. 
This result seems particularly useful since the common lore was that the presence of magnetic fields 
induce a singular modes. As discussed, this presumption was a simple consequence of the fact 
that, in most studies, the anisotropic stress of the neutrinos was simply neglected. In the present 
analysis, on the contrary, the four-fiuid system describing the most general initial conditions for 
CMB anisotropies has been consistently studied and solved. 
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The present paper opens various possible developments. First of all the analysis reported in 
the present paper involves the lowest order in the tight coupling expansion. It will be interesting 
to go beyond and study, within the same framework, the effects arising to first-order in the tight 
coupling expansion. This will allow the treatment of various phenomena (including Faraday effect) 
that have been discussed, up to now, just assuming the standard adiabatic mode that is, however, 
itself modified by the presence of the magnetic field as demonstrated in the present paper. Numeri- 
cal studies on the relevance of magnetized fluctuations correlated (or anticorrelated) with adiabatic 
and non-adiabatic modes are now conceptually possible since the super-horizon curvature pertur- 
bations have been carefully computed (after equality) starting from the solution of the system of 
perturbations after neutrino decoupling. Possible strategies of parameter extraction adapted to 
this new situation should be discussed in future investigations. 



35 



A Power spectra 

In the present appendix the main conventions for the evaluation of the power spectra wih be 
summarized. This technical discussion is relevant since there are conventions that are used in the 
magnetic field literature (for instance concerning the relative normalization of the power spectra) 
and that differ from the conventions used in CMB physics. According to the conventions adopted 
in this paper, the Fourier components of the curvature and entropy fluctuations are defined as 

nx. r) = j^^ J d^kn{k, r)e-*^^-; S{x, r) = -^^ J d^kS{k, r)^^'^ , (A.l) 

where TZ*{k,T) = TZ{—k,T) and 5*(fc,r) = 5(— fc, r). The correlation functions of the primordial 
adiabatic and isocurvature modes are then defined as 

{n,ik,T)n,ip,T)) = ^Vnik,T)6^^Hk + p), {S,ik,T)SUp,T)) = ^Vsik,T)6^'Hk + p}, 

(A.2) 
where k = \k\ and where the power spectra can be parametrized as 

^ \ rir-l / k ^ "i-l 



rn{k,T)=An[^-j , rsik,T) = As[^-j . (A.3) 

In Eq. (A.3), kp is the "pivot" scale at which the spectra are normalized. A typical choicc"^^ for 
kp is 0.05 Mpc^"*^. If we allow a correlation between adiabatic and isocurvature modes, the initial 
value problem may require a further primordial correlation function, namely: 

2^2 



{n,{k,T)SMr)) = -^\/'Pn{k,T) Vs{k,T) cos jrsS^^Hk + p), (A.4) 

where jrs is the correlation angle. If cos7r.s = the primordial modes are not correlated. 
With similar notations, the Fourier modes for the magnetic field intensity are 

B^{x,T) = ^-l^Jd^kB,{k,T)e-'^'^-^, (A.5) 

again with B^{k,T) = Bi{—k,T). The magnetic fields studied in the present paper are fully inho- 
mogeneous and characterized by their two-point function, i.e. 

27^2 
{B,{k,T)B\p,T)) = —P^{k)Ps{k,T)5^'Hk+P), (A.6) 

where 

Pnk)=(si-^~^, PBik,r)=A^^^y'\ (A.7) 



^^ Another possible choice for kp is 0.002 Mpc ^. See [3] and also [8] for a discussion of the sensitivity of the 
parameter estimation on kp. 
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The amplitude of the magnetic power spectrum, ^B) can be traded for the magnetic energy density 
regularized over a typical scale of size L. In fact, by using a Gaussian window function e~ " for 
each magnetic field intensity, the (real space) magnetic autocorrelation function at two coincident 
spatial points is 

i?i = (S,(x,r)SXx,r))=AB(^)"'V(^) (A.8) 

where ki = 27r/L. Equation (A.8) implies that the amplitude of the magnetic power spectrum 
appearing in Eqs. (A. 6) and (A. 7) can be expressed in terms of B\ as 

•*= = (^'''f|2)(fe)'- <*■'" 

where, for convenience, the magnetic spectral tilt e = {m — 1) has been introduced. Different 
authors (see, for instance, [11, 21, 22, 23]) employ slightly different conventions both for Eqs. 
(A.7)-(A.9). In particular, sometimes the factor l-ii^ jk^ appearing in Eq. (A. 6) is dropped and 
reabsorbed in the definition of PBik) which is, in turn, often parametrized as Pb = Ak"^ (without 
any typical pivot scale). If these two choices are consistently adopted, a nearly scale invariant 
magnetic energy spectrum (see below) occurs for —3 < m < —2. Within the present conventions, 
that are homogeneous with the ones customarily adopted in the analysis of primordial adiabatic 
and isocurvature modes, the nearly scale-invariant magnetic energy spectrum is obtained in the 
range < e < 1. 

The nearly scale-invariant spectrum is particularly interesting for the purposes of the present 
investigation. In this limit, in fact, the magnetic field energy spectrum (as well as the spectrum of 
the Lorentz force) is practically unaffected by diffusive damping provided by the effective Alfven 
damping scale [20, 21]. The Alfven damping occcurs for typical wavelengths that are of the order of 
La = VcsLs where Ves is the effctive Alfven velocity and Ls is the Silk damping scale. Since Ves will 
be proportional, in the present notations, to B^/p^, we can also appreciate that La will simply be 
a fraction of the Silk scale. In the case covered by the present investigation, on the contrary, we are 
not interested in setting bounds on very steep (violet) spectra of the magnetic energy density that 
are, in any case, already strongly constrained by the analysis of the vector and tensor modes. On 
the contrary, since in the present analysis, the magnetic fields induce some correlated isocurvature 
modes, it is also reasonable to deepen the analysis in the nearly scale-invariant case. With these 
specifications, the typical scale where diffusion becomes important will not be directly related with 
the Alfven damping but rather with the typical scale at which the wavelength of the fluctuations 
becomes comparable with the photon mean free path corresponding to a comoving wave-number 

Prom Eq. (A. 6) it is possible to compute, by convolution, the appropriate quadratic observables 
appearing in the treatment of the scalar problem. In particular, as discussed in the bulk of the 
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paper, the effect of the magnetic fields on scalar modes depends upon the magnetic energy and 
pressure densities, (i.e. 6pB and 6pb) as well as on the divergence of the Lorentz force. This last 
quantity, in MHD, is proportional to V • [(V x B) x B]. Recalling, in fact, the notations established 
in Eqs. (3.12) and (3.13) we have that, in real space 

^^(^^'") = ^tI^ = (^/^''^ ^B(g»e-^«^ 6ps = ^, (A.ll) 

where 

In Eq. (A.12) the quantity p = p^{T)a'^{T) has been introduced for practical reasons. Consistently 
with Eq. (A. 5), the correlation function of the magnetic energy density is then: 

(OB(fc,r)J^B(p,r)) = '^Vn{k,T)5'^^\k + p). (A.13) 

where 

Vn{k)=He)^lL[YS, (A.14) 

(6-e)(27r)2= ,^ , 

^(-) = .(3-2e)r^(;/2) ' ^^-''^ 

nBL = ^, PBL = ^. (A.16) 

P-y Svr 

Notice that the obtained expression for Vn follows from the mode-coupling integral 

^B(g) = ^^ /d34l + Mf±^U(g + fc)PB(fc). (A.17) 

4(27r)^ p^a"^ J [ k^{q + k)^ J 

in the limit < e < 0.3. According to Eq. (3.12), the correlator of (Tb determines, together with 
Ob the Lorentz force. As far as the two-point function of ctb is concerned, we can start by noticing 
that 



"^^"^ ^ (2vr)3/2 



1^ I d'qe-^'^-^aBiq), (A.18) 

where 

U(g - p)B\p) + 3{q^ - p^)fBj{p}Bi{q- p) 



"^(^) = (2vr)3/n6vrlv. ? ^ ^'"^ 



(A.19) 



Consequently, 

(aB(QlaB(g-')) = ^P.(g)5(3)(g-+ g'). (A.20) 



where, for e < 0.2, 



^,, (3-2e)(188-4£2-66e)^ , , ,,„^, 
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In the general phenonienological situation, the magnetic field may also be correlated with the 
primordial adiabatic and isocurvature modes. In this case, two new cross-correlations have to be 
taken into account when computing, for instance, the Sachs- Wolfe plateau: 



{S4k)nB{p)) = ^\/rs{k) Pn{k)cosTshS^''Hk + p}. (A.22) 



Again, in Eq. (A.22) the various "angles" parametrize the correlation between different modes. If 
Irs = 7rb = 7bs = 7i"/2 the modes are not correlated. 



39 



References 

[1] D. N. Spergel et al. [WMAP Collaboration], Astrophys. J. Suppl. 148, 175 (2003). 

[2] H. V. Peiris et al, Astrophys. J. Suppl. 148, 213 (2003). 

[3] D. Spergel et al. [WMAP Collaboration], astro-ph/0603449. 

[4] J. M. Bardeen, Phys. Rev. D 22, 1882 (1980). 

[5] H. Kodama and M. Sasaki, Prog. Theor. Phys. Suppl. 78, 1 (1984). 

[6] K. Enqvist, H. Kurki-Suonio and J. Valiviita, Phys. Rev. D 62, 103003 (2000). 

[7] J. Valiviita and V. Muhonen, Phys. Rev. Lett. 91, 131302 (2003). 

[8] H. Kurki-Suonio, V. Muhonen and J. Valiviita, Phys. Rev. D 71, 063005 (2005). 

[9] K. Subramanian, arXiv:astro-ph/0601570. 

[10] M. Giovannini, Class. Quant. Grav. 23, Rl (2006). 

[11] A. Mack, T. Kahniashvih and A. Kosowsky, Phys. Rev. D 65, 123004 (2002). 

[12] A. Lewis, Phys. Rev. D 70, 043011 (2004). 

[13] T. Kahniashvih and B. Ratra, Phys. Rev. D 71, 103006 (2005). 

[14] A. Kosowsky and A. Loeb, Astrophys. J. 469, 1 (1996). 

[15] M. Giovannini, Phys. Rev. D 56, 3198 (1997). 

[16] C. Scoccola, D. Harari and S. Mollerach, Phys. Rev. D 70, 063003 (2004). 

[17] M. Giovannini, Phys. Rev. D 71, 021301 (2005); 

[18] A. Kosowsky, T. Kahniashvih, G. Lavrelashvih and B. Ratra, Phys. Rev. D 71, 043006 (2005). 

[19] L. Campanelh, A. D. Dolgov, M. Giannotti and F. L. ViUante, Astrophys. J. 616, 1 (2004). 

[20] K. Subramanian and J. D. Barrow, Phys. Rev. D 58, 083502 (1998). 

[21] K. Subramanian and J. D. Barrow, Phys. Rev. Lett. 81, 3575 (1998). 

[22] K. Subramanian and J. D. Barrow, Mon. Not. Roy Astron. Soc. 335, L57 (2002). 

[23] K. Subramanian, T. R. Seshadri and J. D. Barrow, Mon. Not. Roy. Astron. Soc. 344, L31 
(2003). 

40 



[24] A. Challinor and A. Lasenby, Phys. Rev. D 58, 023001 (1998). 

[25] M. Bucher, K. Moodley and N. Turok, Phys. Rev. D 62, 083508 (2000). 

[26] C. G. Tsagas and J. D. Barrow, Class. Quant. Grav. 15, 3523 (1998). 

[27] C. Tsagas and R. Maartens, Phys. Rev. D 61, 083519 (2000). 

[28] M. Giovannini, Phys. Rev. D 70, 123507 (2004). 

[29] D. Biskamp, Non-linear Magnetohydrodynamics, (Cambridge University Press, Cambridge, 
UK, 1994). 

[30] T. J. M Boyd and J. J. Sanderson, The Physics of Plasmas, (Cambridge University Press, 
Cambridge, UK, 2003). 

[31] N. A. Krall and A. W. Trivelpiece, Principles of Plasma Physics, (San Francisco Press, San 
Francisco, USA, 1986). 

[32] R. Beck, A. Brandenburg, D. Moss, A. Shukurov and D. SokolofF, Ann. Rev. Astron. Astrophys. 
34, 155 (1996). 

[33] F. Govoni and L. Feretti, Int. J. Mod. Phys. D 13, 1549 (2004). 

[34] R. Cassano, G. Brunetti and G. Setti, e-print astro-ph/0601702. 

[35] Y. Xu, P. P. Kronberg, S. Habib and Q. W. Dufton, Astrophys. J. 637, 19 (2006). 

[36] A. Dar and A. De Rujula, Phys. Rev. D 72, 123002 (2005). 

[37] M. Giovannini, Int. J. Mod. Phys. D 13, 391 (2004). 

[38] A. Brandenburg and K. Subramanian, Phys. Rept. 417, 1 (2005). 

[39] A. Lazarian, E. Vishniac and J. Cho, Lect. Notes Phys. 614, 376 (2003). 

[40] J. M. Bardeen, P. J. Steinhardt, and M. Turner, Phys. Rev. D 28, 679 (1983). 

[41] J. M. Bardeen Lectures given at the 2nd Guo Shou-jing Summer School on Particle Physics 
and Cosmology (Nanjing, July, 188) DOE/ER/40423-01-C8. 

[42] J. Hwang, Astrophys. J. 375, 443 (1990). 

[43] M. Giovannini, Phys. Lett. B 622, 349 (2005). 

[44] R. Brandenberger, R. Kahn, and W. Press, Phys. Rev. D 28, 1809 (1983). 

[45] D. H. Lyth, Phys. Rev. D 31, 1792 (1985). 

41 



[46] C. P. Ma and E. Bertschingcr, Astrophys. J. 455, 7 (1995). 

[47] I. Brown and R. Crittenden, Phys. Rev. D 72, 063002 (2005). 

[48] A. G. Riess et al. [Supernova Search Team Collaboration], Astrophys. J. 607, 665 (2004). 

[49] H. Jorgensen, E. Kotok, P. Naselsky, and I Novikov, Astron. Astrophys. 294, 639 (1995). 

[50] W. Hu and N. Sugiyama, Astrophys. J. 471, 542 (1996). 

[51] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions (Dover, New York, 
1972). 

[52] A. Erdelyi, W. Magnus, F. Obehettinger, and F. Tricomi, Higher Trascendental Functions (Mc 
Craw-Hill, New York, 1953). 



42 



